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Abstract 

Recently the ATLAS and CMS collaborations presented preliminary results of 
Standard Model Higgs searches and reported excesses of events for a Higgs boson 
at 124 — 126 GeV. Such a Higgs mass can be naturally realized, simultaneously 
explaining the muon g — 2 anomaly, in gauge-mediated SUSY breaking models 
with extra vector-like matters. Upper bounds are obtained on the gluino mass, 
mg < 1.2(1.8) TeV, and on the extra vector-like quark mass, Mqi < 1.0(1.8) GeV, 
in the parameter region where the Higgs boson mass is 124 — 126 GeV and the 
muon 5 — 2 is consistent with the experimental value at the la {2a) level. The LHC 
prospects are explored in the parameter region. It is found that some of the regions 
are already excluded by the LHC, and most of the parameter space is expected to 
be covered at y/s = 14 TeV. A study on the extra vector- like quarks, especially 
current bounds on their masses and prospects for future searches, is also included. 



1 Introduction 



Recently ATLAS and CMS collaboration have presented preliminary results of Standard 
Model (SM) Higgs searches, in a dataset corresponding to an integrated luminosity of up to 
4.7- 4.9 fb"^ collected at y/s = 7TeV in the LHC [I]. Interestingly, excesses of events are 
observed for a Higgs boson mass hypothesis close to 124 — 126 GeV at both experiments. 
Such a Higgs boson mass is consistent with the prediction of the supersymmetry (SUSY), 
which is one of the best candidates for new physics beyond the Standard Model. 

On the other hand, the supersymmetric standard models are one of the most natural 
scenarios which can explain the anomaly of the muon anomalous magnetic moment (the 
muon g — 2). Latest studies have reported the discrepancy of the measured muon — 2 [2] 
from the Standard Model prediction by more than the 3a level [31 Hj. 

Recently we have shown in Ref. [3] that a relatively heavy Higgs boson and the 
discrepancy of the muon g — 2 can be simultaneously explained in SUSY models with 
vector-like matters, in the frameworks of gauge-mediated SUSY breaking (GMSB) models 
and minimal supergravity (mSUGRA) models. Interestingly, it was shown that the Higgs 
boson mass of 125 GeV can be consistent with muon g — 2 within la in GMSB, and 
within 1.2cr in mSUGRA models]^ The key point is that the vector-like matters coupled 
to the Higgs field can enhance the Higgs boson mass [HI El E] in the same mechanism as 
the top (s)quark [9l [10]. There has been recently growing interest in SUSY models with 
vector-like matters to enhance the Higgs boson mass [5l HU [121 113] • 

In this paper, we extend our previous work on the GMSB models with extra matters in 
Ref. [5] in light of the current Higgs search results. Among various high-energy models, the 
GMSB [H] is one of the most attractive models from the phenomenological viewpoints, 
since dangerous flavor-changing processes and CP violations are naturally suppressed. 
However, the simplest GMSB models face the difficulty of explaining the mass of the 
Higgs boson in the range of 124 — 126 GeV, since the scalar trilinear coupling of the top 
quark is small and its contribution to the Higgs potential is suppressed. In contrast, the 
GMSB models with extra vector-like matters [5l |12] are one of the most attractive and 

^Here and hereafter, the "Higgs boson" refers to the lightest CP-even Higgs boson in the minimal 
SUSY Standard Model (MSSM), which is the Standard- Model-like Higgs boson. 
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phenomenologically viable SUSY models which can realize the Higgs mass 124 — 126 GeV. 
In this paper, prospects of SUSY discovery in the LHC are explored in light of the current 
Higgs search results and the muon g — 2 discrepancy. In the region where the Higgs boson 
mass is 124 — 126 GeV and the muon (7 — 2 is consistent with the experimental value at 
the Icr (2cr) level, the gluino mass is bounded as < 1.2(1.8) TeV. The LHC signature 
depends on species of the next-to-lightest SUSY particle (NLSP), which is either stau or 
neutralino. Both scenarios are studied and it will be shown that the most of the parameter 
region is expected to be covered by the LHC in future at the y/s = 14 TeV collision. 
Furthermore, the masses of the extra vector-like quarks are also bounded from above as 
Mqi < 1.0(1.8) TeV. The LHC search for extra vector-like quarks is also investigated. 

This paper is organized as follows. The model is introduced in Sec. [2j The mass 
bounds on the SUSY particles and the vector-like quark are obtained from the present 
situation of the Higgs boson searches and the muon g — 2 in Sec. [3j In Sec. |4] we discuss the 
LHC prospects. The cases of stau NLSP and neutralino NLSP are separately studied in 



Sec. 4.2 and Sec. 4.3, respectively. Phenomenology of extra vector-like quarks, especially 



current experimental bounds on their mass and further LHC search, is studied in Sec. AA 
Sec. |5] is devoted to summary and discussions. 

2 Model 

We consider the simplest GMSB models, which are parametrized by the messenger scale 
Mjness, the soft mass scale, A = -Fmess/^mess, the messenger number N^, the ratio of the 
Higgs vacuum expectation values tan (3 = (Hy) / {H^) and the sign of the Higgsino mass 
sgn(/i). In addition, a vector-like pair of complete SU(5) multiplets, 10 = {Q\ U\E') and 
10 = (Q', U',E'), is introduced, which has a superpotential, 

W = Y'Q'H^U' + V'Q'HdJ' + Mq,Q'Q' + Mu'U'lJ' + Me'E'E', (1) 

and corresponding soft SUSY breaking terms. We follow Ref. [S] for the definitions and 
conventions. In the following, we set = 1 in order to preserve the perturbativity of the 
gauge coupling constants up to the GUT scale, and sgn(/i) = 1. We also assume Y" ~ 0, 
since the extra down-type quark decreases the SUSY contribution to the Higgs mass if 
Y" is sizable (see Ref. |i5j). 
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The soft SUSY breaking parameters are set by the messenger fields at the messenger 
scale, which are developed down to the weak scale following the renormalization group 
(RG) equations. In the numerical analysis, the RG equations are solved at the two-loop 
level by the SuSpect package [15] which is modified to introduce the vector-like matter]^ 
The Higgs boson mass is evaluated at the NLO level for the minimal SUSY standard 
model (MSSM) contribution by the FeynHiggs package [16], and the contribution from 
the extra matter is included at the one-loop level. The estimation of the muon g — 2 is 
obtained from FeynHiggs. See Ref. [S] for more details. 

3 Higgs mass and muon g — 2 

In this section, the Higgs boson mass and the muon g — 2 are studied in the GMSB model 
with the vector-like matter, presented in the previous section. The results are summarized 
in Figs. [T| |2| and [3j The Higgs boson mass and the muon g — 2 are shown as functions 
of the messenger scale and the soft mass scale in Fig. [l| the gluino mass and tan /3 in 
Fig. [2| and the gluino mass and the vector-like quark mass in Fig. |3} It is found that the 
mass of the Higgs boson at 124 — 126 GeV is easily realized with the muon g — 2 anomaly 
explained. Importantly, we obtain upper bounds both on the soft parameters and the 
extra vector-like quark mass. In the following, we discuss these features in more detail. 

An extra contribution to the Higgs potential arises due to the extra up-type quark 
which couples to the up-type Higgs field. It becomes significant when the Yukawa coupling 
of the Q'-Hu-U' interaction is as large as ~ 1. This is guaranteed by the RG evolution, 
since an infrared fixed-point exists at F' ^ 1 [8J. The contribution to the Higgs potential 
is similar to that of the top (s) quark and is enhanced when there is a large hierarchy 
between the scalar and fermion masses of the extra up-type quark. Although a large 
trilinear coupling of the extra up-type squark can also enhance the contribution, this 
is quite unlikely to happen because the coupling has an infrared fixed-point at a small 
value IS]. This feature is different from the so-called m/i-max scenario of the MSSM, 
where the trilinear coupling of the stop is tuned to maximize the stop contribution to 
the Higgs boson mass. In summary, the Higgs mass of 124 — 126 GeV is saturated in the 
^The relevant RG equations are summarized in Ref. [5]. 



4 



Mq, = 0.6 TeV Mo' = ITeV 



12 



^11 
> 

O 10 



a 9 



O 



Q' 



Mq' =0.6 TeV Mq, = 1 TeV 





O.O 

\ 


1.0 


1.2 


1.4 


u 






\ 

\ 

\ 

\ 

\ 












Si 

h-' 


\ 


. V 


• 








" A 




\ 

\ 

\ 

\ 










h- ' 




\ 

\ 










O 






\ 








o 






\ 

\ 








O 






\ 

\ 

\ 








ft) 






\ 








< 






\ 

\ 














\ 

\ 

1 » 










40 60 80 100 120 140 160 180 200 40 60 80 100 120 140 160 180 200 

A (TeV) A (TeV) 

Figure 1: The regions of the Higgs boson mass of 124GeV < rrih < 126GeV and the muon 
g — 2 are shown for tan /3 = 30. The two green bands correspond to the Higgs mass for the 
vector-hke quark mass, Mq/ = Mjji = 600 GeV and 1000 GeV, while the region favored 
by the muon g — 2 anomaly is displayed by the orange (yellow) region for the la{2a) level. 
The gluino (squark) masses are also given by the black solid curves in units of TeV in the 
left (right) panel. On the blue dashed line, the mass of the lightest neutralino is equal to 
that of the lighter stau. The stau mass < 100 GeV has been excluded by LEP [T7] if it is 
long-lived. 



following two parameter regions (see Figs. [T}|3]): 

• the soft SUSY-breaking mass scale is large, 

• the SUSY-invariant mass of the vector-like matter is small. 

The SM prediction of the muon g — 2 is currently 3 — 4(T smaller than the experimental 
value. The SUSY contributions to the muon g — 2 is enhanced when the soft SUSY- 
breaking mass scale is low and tan /3 is large. This provides an upper bound on the soft 
masses, which is essential for the study of the discovery potential of the SUSY particles 
at the LHC. In detail, the main contribution comes from the chargino-muon sneutrino 
diagram. In addition, the neutralino-smuon diagram can give a comparable contribution 
when the fi term is large. This is indeed the case in the vector-like matter models, which is 
one of the characteristic features of the models. The feature is because the extra up-type 
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Figure 2: The regions/lines are the same as Fig. [l] except for the black hnes with various 
messenger scales. The black lines denote contours of the production cross sections of the 
SUSY event via the EW processes in units of fb for the LHC at ^/s = 7 TeV. 
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Figure 3: The regions/lines are the same as Fig. [Tjexcept for the black hnes with tan (3 = 20 
(left) and 40 (right). In the right panel, contours of the lightest stau mass are shown by 
the black solid lines in units of GeV. 



quark draws down the soft mass of the up-type Higgs, mj^^, to a large negative value 
during RG running in a similar way as the top quark does. 

In Fig. [l] the Higgs mass and the muon g — 2 are evaluated in the A-Mmess plane 
for tan /3 = 30. The gluino and squark masses are also shown in the figure. The Higgs 
mass of 124 GeV < rrih < 126 GeV is realized in the green bands, where the SUSY- 
invariant vector-like quark masses are Mqi = Mur{= Me') = 600 GeV and 1 TeV. As the 
SUSY-invariant mass increases, the bands shift rightwards because larger soft masses are 
required. The regions which are consistent with the experimental value of the muon g — 2 
at the Icr {2a) level are displayed by the orange (yellow) regions. 

It is emphasized that the Higgs boson mass of 124 — 126 GeV is naturally realized in 
the region where the muon g — 2 is even within the la level [5] . As expected, the muon 
g — 2 prefers a small soft mass scale A. It is also found that given the Higgs boson mass, 
a smaller messenger scale is more consistent with the experimental value of the muon 
g ~ 2. This is because as the messenger scale increases the EW gaugino masses enhance 
the slepton masses during RG evolutions. On the other hand, as noticed from Fig. [ij the 
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gluino mass is less sensitive to Mmess, because the /3-function of the gluino mass vanishes 
at the one-loop level for a pair of extra 10 + 10. On the other hand, a light extra up-type 
quark is also favored, since masses of the extra matter components become hierarchical. 

There are two distinct regions from the viewpoints of the LHC discovery. When the 
lightest neutralino is the NLSP, the LHC signature looks like that of mSUGRA unless 
the gravitino mass is very light. If the stau is the NLSP, it becomes long-lived without 
leaving the missing energy in the detectors as long as the gravitino is relatively heavy. In 
Fig. [TJ the blue dashed line is drawn, where the stau has the same mass as the neutralino, 
m^o = mf^. The stau is the NLSP in the left region of the line, while the neutralino 
becomes lighter rightwards. The LHC prospects will be studies in the next section. 

Contours of the Higgs mass and the muon g — 2 are shown in Fig. [2] as a function of the 
(physical) gluino mass and tan /3 for several messenger scales. As tan /3 is larger, heavier 
SUSY particles are allowed by the muon g — 2. On the other hand, the stau becomes the 
NLSP above the blue dashed line, because the r Yukawa couphng is enhanced by a large 



tan (3. As we will discuss in Sec. |4.2[ the region with stau NLSP is almost excluded by the 
current LHC search. On the other hand, in the region with the neutralino NLSP, the soft 
mass scale, i.e. the gluino mass, is severely constrained by the muon g — 2. For a Higgs 
boson mass 124-126 GeV, the gluino mass is required to be less than 1.2(1.8) TeV for the 
muon g — 2 at l(T(2cr), for Mmess = 10^ GeV. Note that a larger messenger scale tightens 
the soft mass scale more severely. 

The Higgs boson mass is also sensitive to the vector-like quark mass. In Fig. [3} contours 
of the Higgs mass and the muon g — 2 are drawn in the plane of the (physical) gluino mass 
and the vector-like quark mass for a couple of tan /3 with M^ess = 10^ GeV fixed. Upper 
bounds on the vector-like quark mass are obtained from the Higgs mass and the muon 
g — 2. When tan /3 is small, e.g., tan /3 = 20 as in the left panel, the neutralino is the NLSP 
in the viable region]^ and the LHC signature resembles that of mSUGRA. The possible 
maximal value of the SUSY-invariant mass of the vector-like quark is 800(1150) GeV for 
the Higgs boson mass of 124 — 126 GeV and the muon g — 2 within the la{2a) level. 
For a larger tan/3, a larger soft mass (gluino mass) is allowed (see Fig. [2]), and therefore 
a heavier vector quark can be consistent with 124-126 Higgs mass. The maximal value 
^The stau NLSP region in the left panel of Fig. [sjis excluded. See Sec. 



4.2 



8 



of Mq/^jji is obtained when the gluino mass takes its maximal value 1.2 TeV (1.8 TeV), 
which leads to Mq,^u' < 1-0 TeV (1.8 TeV). 

On the other hand, the mass spectrum is different in the right panel of Fig. |3| where 
tan /3 = 40. The stau is lighter than the neutralino when the soft mass scale, i.e. the 
gluino mass, is small, whereas the mass relation becomes inverted as the scale increases, 
because the left-right mixing of the stau mass matrix is suppressed compared to the 
chirality-conserving components. From Fig. [2| it is found that the inversion happens at 
around rrig = 1.9 TeV for tan /3 = 40 and M^^ss = 10^ GeV. Thus, the stau is the NLSP in 
the whole region favored by the Higgs search results and the muon g — 2. We will discuss 
the LHC discovery/exclusion in the next section. 

4 LHC prospects 

In this section, the LHC prospects for the GMSB models with vector-like matter are 



explored. In Sec. |4.1[ a typical mass spectrum of the model is briefly summarized, paying 
attention to the LHC search. The cases of stau NLSP and neutralino NLSP are discussed 
in Sec. 4^ and 4.3[ respectively. The LHC searches for vector-like quarks are investigated 



in Sec. 1131 

4.1 Mass spectrum 

As discussed in the previous section, the soft mass scale is bounded by the muon g — 2, 
and the vector-like quark mass is limited by combining the Higgs boson mass and the 
muon g — 2. The LHC reach for SUSY particles very much depends their mass spectrum. 

Let us first discuss the gravitino mass. In GMSB, the gravitino is the lightest SUSY 
particle (LSP), and the NLSP decays to the gravitino. The LHC signature depends on 
the gravitino mass m3/2 as well as species of the NLSP. The NLSP decay length is given 
by CT ^ (9(10 m)(mNLSp/100GeV)~^(m3/2/lkeV)^. In the region where the Higgs boson 
mass of 124 — 126 GeV and the muon g — 2 discrepancy are simultaneously explained, the 
F-term of the messenger sector satisfies F^^ss = ^messA > (9(10^^GeV^). On the other 
hand, the total SUSY breaking scale, -Ftotai, is typically much larger than F^^ss except for 
the direct GMSB models. Thus, in the following analysis, we assume that the gravitino 



9 



mass, m3/2 — Ftotai/^^p, satisfies 7713/2 ^ 0{keV), and hence the NLSP behaves as a 
stable particle at the LHC. Note also that the gravitino mass of (9(0.1 — l)keV is strongly 
disfavored from cosmological points of view fl8\ [19] . Signatures of in-fiight decays of the 
NLSP will be briefiy discussed in Sec. [5j 

The SUSY particles relevant for the discovery depend on the signatures of the SUSY 
events. When the neutralino is the NLSP, jets with a large transverse momentum as well 
as a large missing energy are responsible for discriminating the SUSY events from the 
SM background. Then, the relevant SUSY channels are productions of the colored SUSY 
particles. Their soft masses are mainly controlled by the gluino mass, M3. When a pair 
of the 10 + 10 multiplet is introduced, the /3-function of the SU(3) coupling constant 
vanishes at the one-loop level and stays large between the weak and messenger scales. 
Thus, the /3-function of the gluino mass starts from the two-loop level. It turns out that 
the gluino mass decreases by ~ 40% by the RG running if the mediation scale is the GUT 
scale. Since the gluino mass is not small at the messenger scale compared to the value at 
the weak scale, the soft masses of the squarks receive more contributions from M3 during 
RG and are raised at the weak scale compared to the case of the MSSM. In the analysis, 
the soft parameters are developed at the two-loop level, which is important for the LHC 
study of the SUSY production. 

When the stau is the NLSP, the SUSY signal is very clean against the SM background, 
and all the SUSY production processes become responsible for the discovery. As the soft 
mass scale increases, the production cross sections of the colored SUSY particles drop 
more rapidly than those of the charginos and/or neutralinos. They are composed of the 
Bino, Wino and the Higgsinos. As mentioned in Sec. [3} the Higgsino mass tends to be 
large in the vector-like matter models. It can be checked that the fi parameter becomes 
quantitatively comparable or larger than the gluino mass. Hence, the lightest neutralino is 
almost the Bino, and the next-to-lightest neutralino and the lightest chargino is dominated 
by the Wino component. 

In summary, the typical relation among the soft parameters is 

m^^3, /i > M3 ~ 3.5M2 ~ 7Mi , (2) 
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which leads to 

where {xf) denotes the neutrahnos (charginos). We have checked that this relation 
indeed holds in the parameter regions of our interest. See e.g., Table [T] and [2] 

4.2 Stau NLSP 

When the stau is the NLSP, the SUSY signal is different from the signature of the 
neutralino NLSP. Such a light stau is realized especially when tan/3 is large because 
the r Yukawa coupling is enhanced. When the gravitino is much heavier than (9(keV), 
the stau is long-lived enough to behave as a stable particle and to leave a track in the 
detectors. Since the stau mass is C(IOO) GeV, its velocity has a wide distribution. The 
signal can be distinguished from the SM backgrounds by choosing slowly-propagating 
staus with proper cut conditions. In fact, the dominant background, which is from the 
muon, can be suppressed significantly by selecting the slow staus with a large transverse 
momentum. 

The ATLAS and CMS experiments published results of searches for the heavy long- 
lived charged particle [201 EI]- They basically follow the above strategy to select the events. 
The events are triggered by the muon system. The clean separation of the stau and SM 
backgrounds is achieved by selecting candidates with a large transverse momentum pt and 
a low velocity /3, focusing on a rate of the energy loss through ionization and measuring 
a time-of-fiight. 

The ATLAS has analyzed the data up to 37 pb^^ at ^/s = 7 TeV, providing a 
lower bound on the stau mass of 110 GeV by the SUSY production via the electroweak 
processes in GMSB |20]. Since the integrated luminosity used in the analysis is limited, 
the constraint is not yet restrictive for the current models. 

The CMS analyzed the data of the integrated luminosity 1.09 fb^^ at ^/s = 7 TeV |21j . 
The long-lived stau was searched for in GMSB, and the bound on the cross section is 
obtained as o" < 3(5) fb for rrif = 300(200) GeV. In this region the main production 
channels of the SUSY events are via the EW processes, i.e. the slepton pair and the 
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Figure 4: The stau velocity distribution with the cuts, 40 GeV < ^^(mf) < 1 TeV and 
\r){mf)\ < 1.5 (left), and the transverse momentum distribution with the cuts, 0.4 < 
P{m.f) < 0.9 and \ri{mf)\ < 1.5 (right). The model points, LLPl (red line) and LLP2 
(blue line), and compared with the GMSB point (yellow region). 



chargino/neutralino productions. As a result, the CMS obtained a 95% CL. lower limit 
on the stau mass > 293 GeV, when the cross section is estimated at the LO. 

In order to apply the CMS bound to the vector-like matter models, the distributions 
of P and pt are compared to the GMSB in Fig. |4j Here, based of the geometry of 
the detectors, the cut on the pseudo rapidity is imposed for the stau. The events are 
generated by using PYTHIA6 [22]. The vertical axis denotes the event number, which 
is normalized by the cross section. The relevant mass spectra of the model points are 
summarized in Table [l] The GMSB point is chosen as (A, Mmess, tan /3, sgn(/i), A^s) = 
(97.9 TeV, 160 TeV, 10, +1, 3), which provides the current CMS bound. The points of 
the vector-like matter models, LLPl and LLP2, correspond to the model parameters, 
(A,M^ess,tan/3,sgn(/i)) = (200 TeV, 10^ GeV, 40, +1) and (A, M^ess, tan /3, sgn(/i)) = 
(205 TeV, 10^ GeV, 47.5, +1), respectively, and give the production cross section of ~ 5 fb 
at the LHC of ^/s = 7 TeV. Here, only the EW processes are considered and the cross 
sections are estimated at the LO. Whereas, the productions of the colored SUSY particles 
are irrelevant for the discovery/exclusion in the vector-like matter models, since they are 
heavy. From Fig. |4| it is seen that the /3 and pt distributions are less sensitive to details of 
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Table 1: The mass spectrum of the GMSB + vector-hke matter models, LLPl and LLP2, 
which is compared to the GMSB point providing the current CMS bound. All masses are 
in units of GeV, if not otherwise stated. 

the model except for the mass hierarchy between the stau and the chargino. Although the 
stau velocity distribution tends to be close to /3 = 1 and the selected event number may 
decrease when the mass hierarchy is enhanced, reduction of the event number is expected 
to be within a factor in the parameter region of our interest. Consequently, it is safe to 
read off the constraint simply by estimating the SUSY production cross section. 

The stau production cross section is constrained by CMS as a < 3(5) fb for = 
300(200) GeV. In Fig. |2| contours of the production cross sections of 3, 5 and 10 fb are 
drawn by the solid black lines. It is found that when the stau is the NLSP (above the blue 
dashed line in the figure), the cross section is larger than 3 fb in the parameter region where 
the muon g — 2 is consistent with the experimental value at 2a for tan /3 < 50. Thus, the 
models with the stau NLSP are almost excluded by the search for the long-lived charged 
particle in CMS or will be accessed soon at the LHC. The constraint becomes tighter when 
the messenger scale is higher. Also, since the stau is the NLSP and its production cross 
section is larger than 5 fb in the yellow region in the right panel (tan (3 = 40) of Fig. [3| 
the CMS result is considered to constrain the region already. Full detector simulations 
are needed to conclude the exclusion definitely. 

4.3 Neutralino NLSP 

In this subsection, we discuss the LHC signatures in the case of neutralino NLSP. We 
assume that the decay length of the neutralino NLSP is much longer than the detector 
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Table 2: A comparison of the mass spectrum of the model point PI and the mSUGRA 
point PI'. The model point PI is shown in Figs. [2] and [3j The model point 
PI in GMSB model with vector-like matters is defined by (A, Mmess, tan /3, sgn(/i)) = 
(95TeV,106GeV,20,+l) and Mq,^u',e' = 600 GeV, while the mSUGRA point PI' is 
defined by (mo, M1/2, Aq, tan/3, sgn(/i)) = (1245 GeV, 355 GeV, 0, 20, -hi). AH masses are 
in units of GeV. 

size. In such a case, the neutralino NLSP escapes from the detectors, leaving a missing 
transverse energy. The signature is similar to the case of neutralino LSP scenarios, which 
has been extensively studied particularly in the context of mSUGRA models. 

In order to demonstrate that a typical LHC signature of the present case is very 
similar to that of the mSUGRA models, we compare a model point in our setup 
to that in mSUGRA, where the gluino mass and squark masses are similar to each 
other. Their mass spectra are shown in Table. |2} The model point PI in our 
setup is defined by (A, M^ess, tan/3, sgn(/i)) = (95 TeV, 10^ GeV, 20, -hi) and Mq,^u',e' = 
600GeV, while the mSUGRA point PI' is defined by (mo, M1/2, Ao,tan/3, sgn(/i)) = 
(1245 GeV, 355 GeV, 0, 20, +1). 

The dominant SUSY events (after typical cuts) are produced by gluino pair production, 
or gluino-squark pair production. After produced, they cause cascade decays: 

(q ^) 9 ^ ggX^(orx?^2) ^ ••• (4) 

where the q and q are quarks and anti-quarks, respectively. A typical SUSY event thus 
includes four jets with a large missing transverse momentum carried by the lightest 
neutralino Xiy which is similar to the case of mSUGRA models. 



14 



1000 




500 1000 1500 2000 2500 3000 3500 4000 
meff [GeV] (after other cuts) 



Figure 5: Mes distribution of the model points PI and mSUGRA point PI' (see Tables. |2| 
for 14TeV LHC with an integrated luminosity lfb~^. 



In Figs. [5] and [6} we compare typical event distributions of the two model points at 
14TeV and 8TeV LHC, respectively. We adopt the following cuts (cf. [23]): Q 

• At least one jet with px > 100 GeV and |?7| < 2.5. 

• At least four jets with px > 50 GeV and |?7| < 2.5. 

• Missing transverse energy E!^^^'^ > 100 GeV. 

• Ef^^ > 0.2Meff. 

• A0(jetj - E^^"""") > 0.2 for i = 1,2,3 (three highest pr jets within \r]\ < 2.5). 

• No lepton (e and fi) with pr > 20 GeV and |?7| < 2.5. 

The distributions of Meg after these cuts are shown in Figs. |5] and [6j Here, the effective 

• Note that these cuts and the event distributions are adopted just to illustrate the similarity of the 
LHC signatures between the two model points. The number of events may be below that of the expected 
background, and more optimized cuts and/or more luminosity may be necessary to exclude/discover these 
model points. 
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Figure 6: The same as Fig. [s} but for 8TeV LHC with an integrated luminosity 15 fb ^. 



mass Meff is defined as 

4 

Meff = Er' + ^pj^*'' (5) 

i 

where the sums run over four highest pr jets within |?7| < 2.5. We used the PYTHIA6 to 
study the kinematics, and the gluino and squark production cross sections are estimated 
by Prospino at the NLO level, which leads to a{pp — > gg, qg, qq) ~ 1400 fb at 14 TeV 
LHC for both model points, and a{pp — > gg, qg, qq) ^ 73 fb (76 fb) at 8 TeV LHC for the 
model point PI (PI')) respectively. Contributions of the other SUSY production channels 
(e.g., the chargino production) are negligible in the signal regions. The LHC detector is 
simulated by the PCS package p^ . 

As can be seen in these figures, the event distributions are similar between the model 
points PI and PI', and the total number of events in the signal region are also comparable. 
It is expected that such a similarity between models holds as far as the gluino and squark 
masses are close, since the SUSY production cross section (after cuts) and the energy 
scale of the jets and missing energy are mostly determined by their masses. 

The signals of mSUGRA and the present model may be different in different analysis 
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(e.g., those with less number of jets and/or with leptons). However, we expect that 
approximate LHC reach can be estimated by comparing the model points with similar 
gluino and squarks masses. 

The latest results from 7 TeV LHC with an integrated luminosity of about 1 fb~^ [25| 
put lower bounds on the gluino mass mSUGRA models. At the ATLAS, the gluino 
mass of < 700GeV is excluded for a squark mass < 1200GeV in the four jet channel |25j . 
In the simplified models containing only squarks of first two generations, a gluino and 
a massless neutralino, the gluino mass below ~ 700GeV is excluded. At the CMS, the 
mSUGRA models with the gluino mass of ~ 750 GeV is excluded for the squark mass of 
> 1300 GeV [26] 1^ Therefore, it is expected that the parameter ranges with gluino mass 
of < 700 GeV in Figs, [l] and |2] have been already excluded. 

As discussed in Sec. [3} in our model, in the region where the Higgs boson mass is 
124 — 126 GeV and the muon g — 2 is consistent with the experimental value at the la 
(2(t) level, there is an upper bound on the gluino mass, rrig < 1.2 TeV (m^ < 1.8 TeV). 
When the gluino mass takes its maximal value around 1.2 TeV (1.8 TeV), the squark mass 
is about 2 TeV (3 TeV). In mSUGRA models, such a model point can be reached at the 
14 TeV LHC with an integrated luminosity of 10 fb'^ (100 fb-^). (See, e.g., |27].) Given 
the similarity of the LHC signatures between the two class of models, it is expected that 
the whole parameter space, where the Higgs boson mass is 124 — 126 GeV Higgs boson 
and the muon g — 2 is explained at the la {2a) level, will be covered at the 14 TeV LHC 
with lOfb-i (lOOfb^i). 

4.4 Vector-like quark search 

Lastly we focus on searches for the vector-like quarks. This is quite an interesting topic 
because these particles are not only peculiar to our model but also relatively light and 
thus within the reach of the LHC. 

We will first review the masses and decay modes of the vector-like quarks. Then 
current experimental bounds on the mass of those particles are discussed. After that we 
will mention prospects of further searches. 

^The CMS exclusion is not necessarily from the four-jet cuts. 
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We set Mq' = Mjji and Y" = for simplicity. Y' is set to be Y' = 1.05, the fixed-point 
value. Also rrih = 125 GeV, rrit = 173.1 GeV and rrih = 4.5 GeV are used in the following 
discussion. 

4.4.1 Masses and decay modes 

This model has three massive vector-like quarks. One of them is down-type, b', and the 
others are up-type, t[ and ^2 with m^^ < m^/^. They must be mixed with the SM quarks 
since otherwise the lightest one would be stable, which is cosmologically disfavored. 

We assume that the vector-like quarks are mixed only with third generation quarks 
{Qs, U3, and D2). The superpotential corresponding to the mixing is 

W^mix = euQsHuU' + e'uQ'Hjj^ + SDQ'HaDs. (6) 

The mixing mass terms m'Q'Qz and rnl'UzU' are absorbed to Mqi and Mu' by redefining 
Q' and U' , without loss of generahty. The size of mixing is assumed to be large enough 
to avoid the cosmological problem and also not to be observed as a heavy stable charged 
particle, but so small that the decay of the heavier vector-hke quarks to the SM particles 
are suppressed compared to its decay to a lighter vector-like quark. 
The mass terms in the Lagrangian are 

-CD (g; U' tK) lY'vsinf] 

\e'jjv sin ^ m J \tL J (7) 




where f ~ 174 GeV is the vacuum expectation value of Higgs. The masses of the vector- 
like quarks are 



mt{,t'^ - Mq, ^1 ± a + , rub' ~ Mq/, 



(8) 



where the mass splitting is characterized by 



Y'vsinp 91.4 GeV f Y' \ 
2Mq, Mq, \l.05j ^ ^ 

with a large tan /3. Smaller Mq, = Mu, leads smaller mass separation. 
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0.01 
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Table 3: Benchmark points for the mixing parameters. The shown values as the branching 
ratios are calculated at m^/^ = 400 GeV, Y' = 1.05, tan/3 = 30 and rrih = 125 GeV; they 
have nontrivial dependence on m^/^ , but are almost stable under the absolute values of the 
mixing parameters as long as the mixing parameters are much smaller than (9(1). 

Next let us consider the decay cascade of the vector-like quarks. Because of our 
assumption that the mixing is small, the possible decay channels of the heavier vector- 
like quarks are summarized as 

^2 ^ b'W, t[h, t[Z, b' t[W, t[ bW, th, tZ, (10) 

where some of them may be kinematically forbidden if the mass separation is smaller. 
b' t[W is forbidden if mt' < rrit'^ + m^y, or Mq/(= Mw) < 370 GeV, and there b' 
decays into 3-body. For the decay of t'^, tZ and th channels are respectively closed below 
Mq, ~ 343 GeV and 378 GeV. 

What is important in LHC phenomenology is that the decay branching ratio of the 
lightest vector-like quark, t\, is determined by the mixing parameters e^/, e'u and IS]- In 
order to promote the discussion, we pick up several mixing patterns shown in Table. |3]as 
benchmark points. The branching ratio is determined by the ratio among the parameters, 
and is insensitive to their absolute values as long as they are much smaller than (9(1). 

At the benchmark point (A), t'^ exclusively decays into bW , or Br(t']^ — )■ bW) = 1 for 
any m^'^. However, except for that point, the decay branching ratio of t[ has nontrivial 
dependence on the mass of t[. (cf. Figs. [7]-[9}) This is mainly because the t[ — > tZ 
and t[ — )■ th channels are closed if the mass of t[ is below the thresholds. Especially, if 
m^/ < nit + "tTT-z — 264 GeV, only the t[ — )■ bW channel is open regardless of the mixing 
parameters. 
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Figure 7: The CMS experimental 95% CL upper limit on the t'i' pair-production cross 
section with 1.14 fb~^ data, assuming t'-quark exclusively decays via t' — )■ tZ, as a function 
of rriti (the red solid line) [28]. The black solid line is the NLO total cross section of t'V 
production, or in other words, the cross section with an assumption that Br(t' — )■ tZ) = 1. 
Considering the branching ratio Br(t'^ — > tZ), this limit may give an upper bound on the 
mass of t[ quarks in our model; we show the corresponding t[t[ cross section with the 
branching effect at the benchmark points as dashed and dotted lines. Note that the line 
corresponding to the point (A) is not shown since BT{t[ — )■ tZ) = at the point. 




Figure 8: The same as Fig. [7| but here t'-quark is assumed to decay exclusively via 
t' — )■ bW channel [3ll |32] • The red and blue solid lines denote upper bounds; see the text 
for detail. The line for (A) overlaps the total cross section line since Br{t[ — )■ bW) = 1. 
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4.4.2 Current experimental bounds 

The search for the vector-hke quarks is similar to that for the fourth generation quarks. 
However no bound can be extracted from current experimental results for the heavier 
vector-like quarks, and 6', since they decay into a lighter vector-like quark. On the 
other hand, the mass of t'^ may be limited from searches for the fourth generation t', 
where it is assumed that t' is directly produced in t't' pair and decays exclusively via a 
specific channel]^ 

Assuming the exclusive decay t' — )■ tZ, the CMS collaboration obtained a mass bound 
mt' > 475 GeV at 95% CL with an integrated luminosity of 1.14 fb~^ |28], as shown in 
Fig.|7l Inthefi gure we also show a{pp t[t'i) x [Br(t']^ — > tZ)]"^ for each benchmark point. 
The total cross section a{pp — )■ t'i^t[) is calculated with HATHOR [22] at the NLO level, 
using CTIO [30] parton distribution function. It gives constraints only for the points (D) 
and (E); m^'^ ~ 300 GeV and rrit'^ ~ 265-325 GeV are excluded at 95% CL, respectively. 
The bound for (E) is obtained because Br(t'^ — )■ tZ) is large, while that for (D) is because 
the t'l — )■ th channel is closed and BT{t[ — )■ tZ) is enhanced for mf^ < 298 GeV. Thus, the 
search for t' — )■ tZ has little sensitivity to Point (D). The bound for (E) will be tightened 
by accumulating more data. 

Two results of the search for t' quarks with Br(t' — > bW) = 1 are also presented by the 
CMS collaboration [311 ES]. The result is shown in Fig. |8} The red solid line "di-lep" in 
the figure, comes from the events in which both of two W bosons decay leptonically (ez/g 
or /iz/^), where 1.14 fb~^ data is used [21]. The blue solid line "semi-lep", is extracted from 
the events with semi-leptonic W decay, i.e. one of the W bosons decays leptonically and 
the other decays into jets [32], using 0.573-0.821 fb^^ data. Among the benchmark points, 
the point (A), where Br(t — )■ bW) = 1, receives a 95% CL lower bound m{t[) > 450 GeV. 

The ATLAS collaboration also presented a result of the search for t' quarks decaying 
via t'i' — qWqW — > qluqW with q = u,d,c,s,b using 37 pb^^ data [33]. Because of the 
looser mass bound due to less integrated luminosity, we do not include this bound in 
figures. 

Now let us focus on the Tevatron experiments. Both the CDF and DO collaborations 

^ Here and hereafter, t' denotes a fourth generation up-type quark, while t[ does the up-type vector-hke 
quarks of the present model. 
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Figure 9: The same as Fig. [Sj but for the Tevatron experiments, t'-quark decay is assumed 
to be via t' — )■ hW channel, and the cross section corresponds to 1.96 TeV pp colhsion. 
The red-sohd, red-dotted, and blue-sohd hues denote the upper bounds; see the text for 
detail. The other lines are the same as Figs. [7]-[8| but for Tevatron cross section. Note 
again that the line for (A) sits on that of the total cross section. 

published mass bounds on t' quarks decaying via the t' — )■ hW channel. As shown as the 
red-solid line in Fig. [9| the CDF collaboration [34J gives a bound m^/ > 358 GeV at 95% 
CL. They also presented a bound for t' quarks decaying via t' — )■ q^W channel, where 
is a generic down-type quark in the SM (the red-dotted line in the figure). The blue-solid 
line shows the bound for t' — t- qdW obtained by the DO collaboration [35J- 

In summary, we would like to emphasize that the Tevatron experiments give mass 
bounds for all the benchmark points. This is because t'^ whose mass is less than mt+mz = 
264 GeV decay exclusively via the channel t' — )■ bW as is discussed above, and thus the t[ 
quark below 264 GeV is excluded by the experiments regardless of the mixing parameter, 
or the branching ratio. However, as is clear from the results of the LHC experiments 
especially for Point (D), it is difficult to tighten the mass bound only from the searches 
for f — 7- tZ and t' — )■ bW. Now we are ready to discuss prospects for further search. 
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Figure 10: The production cross sections of pp{pp) — )■ t[t[ as functions of the t'l quark 
mass, calculated with HATHOR [29] at NLO level. 



4.4.3 Prospects of further searches 

First let us review prospects for the search capability at the benchmark points (A)-(E). 
Currently, all of the benchmark points receive lower bounds on the mass from the t' — )■ bW 
search by the CDF experiments as shown in Fig. [9j However, as these bounds owe much to 
the closure of the other decay channels, further search on this channel has little sensitivity 
for points (D) and (E). For the point (E), where t'l — )■ tZ is the dominant decay channel, 
the t' — )■ tZ search will give a bound soon. On the other hard, the point (D) needs 
searches for t' quark with t' — )■ th decay channel, for which no result is published yet. 

As the branching ratio depends largely on the ratio among the small mixing parameter, 
it is of great importance to search for all of the decay modes, especially t[ — )■ th channel. 
This channel is promising for t' search, because the Higgs boson around 125 GeV is 
expected to decay into bb channel and thus three 6-quarks are expected from one t'l, 
or more than three in an event. This characteristic signature would reduce background 
events and, up to 6-tagging efficiency, a good signal over background ratio is expected. 

Search for the heavier vector-like quarks, b' and is also well worth doing. These 
quarks finally decay into one (or three) b quark(s) and at least one vector boson, so there 
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are at least two b quarks and at least two vector bosons. Since the event has so many 
particles, analyses with detector simulation are important for future studies. 

The LHC upgrade to 14TeV is expected to provide good sensitivity to the vector- 



like quarks. As shown in Fig. 10, the t'i' production cross section would be more than 
ten times larger than the 7TeV LHC. This will be delightful especially for the heavier 
vector-like quark search. 



5 Summary and discussion 

We have investigated the GMSB models with vector-like matters, paying particular 
attention to the Higgs mass, muon g — 2, and the LHC discovery prospects of SUSY 
particles. In the region where the Higgs boson mass is 124-126 GeV and the muon g — 2 
is consistent with the experimental value at the la (2cr) level, there is an upper bound 
on the gluino mass, nig < 1.2 TeV {nig < 1.8 TeV), and an upper bound on the extra 
vector-like quarks, 1.0 (1.8) TeV. The NLSP is either stau or neutralino, and the LHC 
prospects of both cases have been studied. Some parameter regions are already excluded 
by 7 TeV LHC, and 14 TeV LHC is expected to cover most of the parameter space. The 
LHC search for extra vector-like quarks was also discussed. 

In this paper, we have mainly investigated the case where the NLSP is long-lived. 
Let us briefly discuss the case where the gravitino is very light, and hence the NLSP 
decays inside the detectors. In the case of the stau NLSP, in-flight decays of the stau 
leave kink signatures. For instance, in the model points LLPl an LLP2 in Table. [T| the 
gluino mass is about 1.8 — 1.9 TeV and the stau mass is about 230 — 280 GeV. According 
to the study in Ref. [36|, for such a mass spectrum, more than 10 kink events can be 
observed at the 14 TeV LHC with an integrated luminosity 10 fb^^, for the stau decay 
length CT ~ 100-4000 mm. For the stau mass 230 — 280 GeV, this corresponds to the 
gravitino mass 0.6 — 6keV. Since LLPl and LLP2 are the points with a relatively large 
gluino mass within the parameter region where a 124 — 126 GeV Higgs boson and a muon 
g — 2 are simultaneously explained, we expect that most of the parameter space can 
be covered at 14 TeV LHC also in the case of in-flight stau NLSP decay. In the case of 
neutralino NLSP with a very light gravitino, typical LHC signatures will be a non-pointing 
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photon [27] and /or a neutralino in-flight decay into Z-boson [3H] • Since there is an upper 
bound on the gluino mass rrig < 1.2 TeV (1.8 TeV) in the case of neutralino NLSP once a 
124 — 126GeV Higgs and the muon g~2 constraints at la {2a) are imposed, it is expected 
that a large part of the parameter space is reached also in this case. 

Concerning the vector-like quark, we again emphasize the importance of the search 
for t' —7- th, the fourth generation quarks decaying into a t quark and a Higgs boson, 
especially if the Higgs boson is really discovered around 125 GeV. The signature of this 
decay channel, three b quarks and a W boson, could be observed in the LHC, and thus 
expected to give a tight bound on the parameters on this model. 

The present model is one of the most attractive and phenomenologically viable SUSY 
models which can explain the Higgs mass 124-126 GeV. Perturbative coupling unification 
is realized. Dangerous fiavor/CP violating soft terms are naturally suppressed, while the 
muon g — 2 can be explained. The model is also cosmologically viable. The gravitino 
can be the dominant component of the dark matter with a reheating temperature Tr ~ 
108GeV(m3/2/lGeV)^i [SQjQln such a case, non-thermal leptogenesis [ID] can explain the 
cosmological baryon asymmetry. Furthermore, the notorious inflaton-induced gravitino 
problem [H], which excludes most of the inflation models in gravity-mediated SUSY 
breaking models, can also be avoided in the present model. 
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